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Abstract 

The dissipative dynamics of a Josephson junction in the Bose-gases is considered within the 
framework of the model of a tunnehng Hamiltonian. The effective action which describes the 
dynamics of the phase difference across the junction is derived using functional integration method. 
The dynamic equation obtained for the phase difference across the junction is analyzed for the finite 
temperatures in the low frequency limit involving the radiation terms. The asymmetric case of the 
Bose-gases with the different order parameters is calculated as well. 

PACS numbers: 03.75.Fi 
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I. INTRODUCTION 



The experimental realization of Bose-Einstein condensation in atomic vapors [1, 2, 3, 4] 
has allowed to observe a great variety of macroscopic quantum effects. In particular, there 
arises a considerable interest to the study of the Josephson effect in the Bose-condensed 
gases as one of intriguing possibilities to explore the macroscopic quantum effects related 
directly to the broken symmetry in the quantum systems. The dynamics of the Josephson 
effect is governed by the difference between the phases of the condensates, playing a role of 
macroscopic quantum variable. 

The theoretical treatment of the Josephson effect includes both the internal effect for 
atoms of a gas in the different hyperfine states and the case of the Bose-condensates spatially 
separated with a potential barrier which acts as a tunneling junction. The latter case due to 
its direct analogy with superconductors seems us more attractive. A lot of work has already 
been done in this direction. In [5] the behavior of the condensate density near the potential 
boundary has been discussed and the quasiclassical expression for the current through a 
potential barrier has been obtained. The articles [6, 7] are devoted to an applicabihty of 
the two-mode approximation in the Josephson junction dynamics. Milburn et al in [8] have 
shown an existence of the self-trapping effect as well as the collapse and revival sequence in 
the relative population. In [9, 10, 11] the nonlinear Josephson dynamics and macroscopic 
fluctuations have been considered, resulting in the optimum conditions [12] to observe the 
Josephson oscillations. Zapata et al [13] have presented a semiclassical description of the 
Josephson junction dynamics. The time-dependent variational analysis of the Josephson 
effect is given in [14] . 

One of the most interesting and important aspects in the Josephson junction dynamics 
from both the theoretical and the experimental viewpoints is the dephasing of the Joseph- 
son oscillations due to coupling between the macroscopic relative phase variable and the 
infinite number of the microscopic degrees of freedom [16, 17]. Historically, in the case of 
the superconducting systems such description of the phase dynamics was developed in the 
middle of 1980's [18, 19, 20]. The most important result was a successive derivation of the 
effective action for the relative phase, revealing the key role of the microscopic degrees of 
freedom in the irreversible dynamics of the superconducting Josephson junctions. From the 
mathematical point of view the response functions in the effective action, which prove to be 
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nonlocal in time, give the full information on the dynamics of a junction. The employment 
of the low frequency expansion for the response functions allows one to obtain the dissipative 
dynamics of a superconducting junction, involving Josephson energy, renormalization of the 
junction capacity (inverse effective mass), and resistance (effective friction) of a junction. 

For the system of two Bose-condensates connected with a weakly coupled junction, it 
is very desirable to trace and explore the dynamics of the relative phase, generalizing the 
method of the derivation of the effective action from the superconducting case to the case 
of the Bose-condensed systems. As we will show in the next sections, the gapless sound- 
like spectrum of low energy excitations in the Bose-condensed gases results in a qualitative 
change of the irreversible phase dynamics compared with that of the superconducting junc- 
tions. So, the main aim of the paper is to derive the effective action for the Bosc point-like 
junction within the framework of the functional integration method in order to find the 
explicit expressions for the response functions and analyze the low frequency dynamics of a 
Bose junction. 

The plan of the article is the following. First, we derive the general expression for the 
effective action depending only on the relative phase for the system of two Bose-condensates 
connected by a point-like junction. Then we consider the case of zero temperature. As a 
next step, we investigate the effect of finite temperature on the phase dynamics. In addition, 
from the low frequency expansion of the response functions we find the Josephson energy, 
renormalization of the effective mass, friction coefficient, and the radiation corrections. The 
latter can be interpreted as a sound emission from the region of a Bose junction. Finally, 
we present the case of an asymmetric junction in the Appendix and summarize the results 
in the Conclusion. 

II. EFFECTIVE ACTION 

First, it may be useful to make some remarks on the geometry of the Bose junction 
and condensates. We keep in mind the case of a point or weakly coupled junction due to 
a large potential barrier between the two macroscopic infinite reservoirs containing Bose- 
condensates. So, we can neglect the feed-back effect of the junction on the Bose-condensates 
and assume that the both condensates are always in the thermal equilibrium state with the 
constant density depending on the temperature alone. The traditional image of such system 
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is two bulks with one common point through which the transmission of particles is only 
possible with some tunneling amplitude. 

So, our starting point is the so-called tunneling Hamiltonian = 1, volume V — 1) 

H^Hi + Hr + H^ + Ht, (1) 

where Hi^j. describes the bulk Bose-gas on the left-hand and right-hand sides, respectively, 

Hi,r = / d'r (- A _ ^ + ^l,r. (2) 

The coupling constant ui^r — ^'nai^r/''^ where as usual a/,r is the scattering length. The 
energy 

H.-'ii^t (3) 



U=[^,+^]e, (4) 



2 V 2 

is analogous to the capacity energy of a junction in the case of superconductors. The 
constant U can be associated with the second derivative of the total energy E = E[Ni, Nj.) 
with respect to the relative change in the number of particles across the junction 

dNf ^ dN^^ 

and usually is estimated as U = {dfii/dNi + dfir/dNr) [13]. In general, it may depend on 
the concrete type of the Bose-junction and simply describes that the energy of the system 
on the whole may depend on the relative number of particles from each bulk. The total 
number of the particles in each bulk is given by 

Ni,, = / d'r (5) 

The term 

Ht^- J dhSr' (r) / (r, r') (r') + /i.e.] (6) 

is responsible for the transitions of particles from the right-hand to the left-hand bulk and 
vice versa. 

To study the properties of the system described by (1), we calculate the partition function 
using the analogy of the superconducting junction 

j VHiV^-^r exp [Se] , (7) 
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/3/2 

/ dr Le, 



(8) 



where the action on the Matsubara (imaginary) time reads 

Se 

Le = JdV {*+ + |: + H. 

To ehminate the quartic term in the action which comes from the "capacity" energy H^, we 
use the Hubbard-Stratonovich procedure by introducing an additional gauge field V{t) on 
the analogy with the so-called plasmon gauge field in metals 



exp 



(9) 



:/Wexp{-/dr 
fVVexp [-JdT 

Next, we follow the Bogoliubov method of separating the field operators into the conden- 
sate and non-condensate fractions, i.e., — Q,r- + ^/,r- Denoting x — (r, r) and introducing 
convenient Nambu spinor notations for the field operators and, correspondingly, matrices 
for the Green functions and tunneling amplitudes, we arrive at the following expression for 
the partition function 



^ VVV^Cexp [So] j p2$exp [S^] . 



(10) 



Thus, we split the initial action into the two parts which correspond to the condensate 

and noncondensate fields, respectively. 



-V) 



2U 



(11) 



S^ = J dxdx' {$+ (G'(o)-i - 7) $ - C+7$ - $+7C} 
For the field operators we used the spinor notations 



V2 



Lr 



, c 




(12) 



In the expression for the condensate part 5*0 we define the amplitude Jq equal to 

7o = y £rd^r' I {r,r'). (13) 

and corresponding to the tunneling process of the condensate-to-condensate particles. It is 
straightforward to obtain the following expressions for the matrix Green functions 



(5(0)- 



(0)-l 











6 {x — x') , 



(14) 



G 



(0)-l 
l,r 



^(0)-l , W,r 
*-^Z,r + ^11 

^02 



77(0)-l ^Z,r I ' 
^Lr + ^11 



(15) 



where the inverse Green functions and self-energy parts are given by the well-known expres- 
sions 



^(o)-i _ d_ _ ^ _ ^ .ZM 



SZ,r Q + spltT + + 

11 — ^'^l,rCi^rCi^rj ^20 — '^l,r Cl^rCl^ri ^02 ~ ''^l,r (^l,r^l,r 

Accordingly, the matrix Green function g\^^ can be represented as 

/ 



(16) 



MO) 



The transfer matrix here has the form 



/ 

r 



I 
7 



Gl,r Fi^r 

, I = I{x,x') = /(r,r') 5 (r-r'). 



(17) 



(18) 



As one can readily see, if we employ a gauge transformation of the field operators 



^i,^ exp [iipi^r (t)] 



(19) 



and impose the conditions (pi — —V/2, ipr = V/2, i.e.. 



(20) 



both normal Gi^r and anomalous Green functions (17) gain additional phase factors with 
respect to the functions in the lack of an external field (notations of [21]). 



Gi,r (r, r') exp {i[ipi^r (t) - (pi,r (r')]) Gi^r {r - r') , 



(21) 



Flt (r, t') exp {i[(pi^r (r) + ifi., (r')]) Fi^r (r - r') . 
The part of action in (11) is quadratic in the non-condensate field operators so we 
can integrate them out. To perform the integration, we employ the well-known formula 



I' 



P^<l>exp -$+q;$ 



= exp 



'p+a-^(3 - Tr [In {a)] 



(22) 
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with a — G^^^ ^ — I and P — IC to arrive at the partition function 

j VipV^Cexp [-S] (23) 

with the effective action given by 



S = So-Tr 



(24) 



In order to run analytically further, it is necessary to make the following approximations. 
First, we expand the second and third terms of (24) in powers of the tunneling amplitude 
/ to the first nonvanishing order. Then, as is stated above, we consider the simplest case 
of a point-like junction putting / (x, x') — IqS (r) 5 (r') 5 {t — t'). The latter also allows us 
to escape from the problem of summing all higher-order terms in the tunneling amplitude 
/, which is inherent in a junction of the plane geometry [17, 22] with the conservation 
of the tangential components of the momentum of a tunneling particle. The problem in 
essence becomes one-dimensional [15, 22] and results in a strongly dissipative low- frequency 
dynamics independent of the tunneling amplitude and governed by the bulk relaxation alone. 
In our consideration this would correspond to the amplitude independent on the x and y 
coordinates, i.e., / oc 5{z)5{z'). Finally, the third approximation we use is a saddle-point 
approximation for the condensate part of the partition function. Substituting ci^r — ^/nol^ 
where no;,r is the density of particles in the condensate fraction, we obtain the expression 
for the partition function depending on the phase difference alone 



Z^= I V(p exp {-Seff [(p]) 



(25) 



The corresponding effective action reads 

Seff[^]^fdT ^ (ff) - 2Ioy/noinor cos (fi 
a{r — t') cos [(/? (r) — (/? (r')] 
+/3 (r - t') cos [p (r) + ip (r')] 



IlSdrdr' < 



(26) 



> . 



Here the response functions can be written using the Green functions 

a (r) = riQrgt (r) + noig+ (r) + G (r) , 

P (r) = norfi (r) + rioifr (r) + F (r) , 



(27) 
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where 



9% (^) = / [Gi, {p, t) ± G^, {P, -r)] , (r) = / ^F,,, (p, r) , 



(28) 



G (r) = 2 (r) ^+ (r) - (r) (r)] , F (r) = 2/, (r) /. (r) . 
The Fourier components of the Green functions of a weakly interacting Bosc-gas arc given 
by [21] 



Gl,r {p,UJn) - 'ti+t^bf' ^^'^ ^^^^^ - ~^1+4M) 



(29) 

(P) = + 2A,,^^p, = £ , \r = 1i/,r-^0/,r- 

Thus, in order to comprehend the dynamics of the relative phase difference across the 
junction, one should analyze the behavior of the response functions a and /? as a function 
of time. 



III. THE RESPONSE FUNCTIONS. 



The calculation of the response functions in the general form is a rather complicated 
problem. However, keeping in view, first of all, the study of the low frequency dynamics of a 
junction, we can restrict our calculations by analyzing the behavior of the response functions 
on the long-time scale. This means that we should find the low frequency decomposition 
of the response functions in the Matsubara frequencies. Next, we will use the procedure of 
analytical continuation in order to derive the dynamic equation which the relative phase </? 
obeys. 



A. Fourier transformation of the a-response. 

From the analysis of the Green function behavior Qi'^i'^) it follows that zero Fourier 
component of the a-rcsponsc function diverges. The simplest way to avoid this obstacle is 
to deal with the difference ci; (a;„) = a (a;„) — q;(0). This corresponds to the substitution 
q;(t) = q;(t) + q;(0)(5(t) into effective action (26) and the second term q;(0)(5(t) yields a 
physically unimportant time-independent contribution into the action, meaning a shift of 
the ground state energy of a junction. Following this procedure, we can arrive at the explicit 
expressions of the Fourier components for all the terms in (27). 
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For the first two terms in the a-response function of (27), we have a simple formula 

- 1 



TTZ/ 

71 



'1 + 



A 



(30) 



and the corresponding expansion to third order in a;/ A is 



{uJn ^ 0) 



TTZ/ 



2^/2 



UJr, 



1 - 



UJr, 



+ 



8 



(31) 



Here, u = mvmA/(\/27r^) is the density of states at the energy equal to A in the normal 
gas. 

In the case of G the calculations of the Fourier components are much more complicated. 
So, we could find the expressions only in the case of zero temperature and the first nonvan- 
ishing order in temperature. Note that nonzero temperature effects are connected with the 

behavior of G{uj) = Go{uj) + Gt{uj). 

For zero-temperature part of G(<^), we obtain 



Gniuj) — —ViVrAlAM I / r . . . TTTT 1 ^ 

' I r I r J J [(A;j;+A^2/)2+a;2](A;j;+A^2/) 



(32) 



^(a:2+l)(2/2+l); • 

Unfortunately, we could not evaluate expression (32) in the explicit analytic form. Thus 
we report here the Fourier expansion up to third order in lo 



Go (c^ ^ 0) 



/A;-A,,\ 1 \uj 



Ai + ArJ 3VAA 



+ 



(33) 



Here we introduced the function 



<^(?) = f + ? (l - - ^1^) - f In 



Li2 
Li2 



1+2 
l-q 



Li2 




+ 



(34) 



where Li2(z) — /°dUn(l — t)/t is the dilogarithm function. The Taylor expansion of (f){q) 
in the case of |g| -C 1 reads: 

Hq)^^-1 + j + --- (35) 
In the case of the equal order parameters A = A; = A^ we obtain: 



Go {cu ^ 0) 



TTt/^A 




2 


\ 5 
.^"2 


1 


u 


+ ... 


8 


A 




~ 3 


A 



(36) 



The other finite-temperature contribution into G can readily be evaluated as 

2 



Gt H 



TTt'/frvA/A^ / 2'kT 



24 




— + w— - 



l Ai + |cj„| 
A, 



' A^ + \ljJn\ 



which in the case of A = A; = A^ gives: 



Gt (co) 



12 



1-1 1 + 



+ ..., 

(37) 

(38) 



Expanding this expression to third order in a;„/A and including zero temperature terms 
yields 



^(c^^^O) 



TTZ/^A 

24 



2 

A 



A 



Svr-f 



1 /27rry 
4 V A J 



(39) 



It is worthwhile to emphasize the appearance a linear term in in Eq.(39) at fi- 
nite temperatures. As we will see below, this results in the manifestation of an additional 

tempcraturc-dcpcndcnt contribution into the dissipation of the junction. The effect can be 
interpreted as a tunneling of normal thermal excitations existing in the system due to finite 
temperatures. 



B. Fourier transformation of the /3-response. 



Like the preceding section, we can evaluate the Fourier components and find the low 

frequency expansion for the anomalous Green functions /(r). For the first two terms in the 
/9-response function (27), we obtain 

3 \ 



3 

+ 8 



^ _5_ 
16 



+ ... . (40) 



/2^1 + |u;„/A| '^"-^0 V2\ 2 

The calculation of F{uin) requires a special attention. In fact, analyzing the formula of the 
Fourier transform 

ttV^i/^Va^a; ^ 1 



F{uJn) 



E 



^ fc=-oo V^(Ai + \uJk\) (A^ +\uJk-UJn\)' 

we see that the zero-temperature expression 

TTiyiiyrV^l^r 7 \ 

J ''^(Aj + |a;'|)(A,+ |a;'-a;|)' 



(41) 



(42) 
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diverges logarithmically at a; = due to behavior of the integrand at large frequencies 
a;' — > oo. To perform integration over cu', we should first pay attention that so far we could 
neglect the dependence on the momentum in the self-energy parts since all the integrals gain 
their values in the region of small momenta. In fact, the self-energy parts depend on the 
momentum and we must use the exact expressions E20 {p) = ^02 (p) = f^oU (p) where U (p) is 
a Fourier component of the interaction between particles, E20 (0) being rioU (0) = iirano/m. 
Here a is the scattering length, uq is the particle density, and m is the mass of a particle. 
The order of the magnitude for the momentum at which U (p) decays can be estimated 
roughly as a reciprocal of the scattering length, i.e., p ~ 1/a. So, we put the upper limit of 
integration equal to the cutoff frequency uOc—l/ ma^ within the logarithmic accuracy. 

Finally, representing F[u!n) as a sum of zcro-tcmpcraturc and finite temperature contri- 
butions F{ijjn) — Fo{u!n) + FT{u!n), wc find the zero-temperature term 



^0 (^n) = 



2 In 



(Va7+a/a7) 



2 In 



■ A, + \/ A/ — |u;„| 



(Va7+Va;) 



+ 



+ arctan 1""!+^'-^'- + arctan \^-\+^^-^^ 



and the finite temperature one, respectively 



Ft M = 



24 

These functions in the case of A; 

FoK) ^TTi/^A I 2 In 



+ 



y'Ar{A;+|a)„|) ^A,(A^+|w„|) _^ 



Ar = A take the form: 



(43) 



(44) 



1 + ^1 + |u;„/A| 



-|- arctan 



|t^n/A| 



2^1 + K/A|^ 



Ft M = 



12 \ IX J \./i + |^„/A| 
Expanding these expressions to third order in a;„/A yields 

/ 

TTI/^A 




(45) 
(46) 



F {Un ^ 0) 



24 



24 In 

2 

A 



+ (¥) - 



3 _ 3 / 27rT Y 
2 i\ A ) 

1 _ 5 ( 27rry 
S \ A ) 



1 



(47) 



Here, as in the case of a-response we emphasize the appearance a linear term in in 
Eq. (47) at finite temperatures which contributes to the dissipation and interpreted as a 
tunneling of normal thermal excitations existing in the system at finite temperatures. 
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C. Functional series of the response functions. 



The Fourier components of the response functions in the form of a series in up to 
third order can be written in the form 



(48) 



Accordingly, the expressions for the response functions in the imaginary time representation 
read as 



a(T) = ai^ {^^Tf)) - «2(^"(t) - as 



2{ttT)* 



3 2 

sii?(7rTT) sin^(7rrr) 



For the sake of brevity, we present expressions for Q;j(T) and /9j(T) in the symmetric case of 
Ai — Ar — A. The general case of A; ^ A^ will be considered in the Appendix. 



+ ..., 



(49) 



sin 



+ .... 



a2 



2 I 1 1 /na£ { 27rry 
2p"'~3V7r { A J 



8A ^ - 



TT 



/9o-7a|i 
Ps 



f-i _ 8 

1^ 3V TT 



12 V A ; 
2 



1 / 27rT\ 

8 V A ; 



(50) 



J. 

2 

37 
8A 



1 + 



57 
16A2 



1 _|_ _2_ / 

15 V TT 



1 -!(¥)■ 

1 _ 5 / 27rTy 

^ 8 V A ; 



Here T/A < 1, 7 = V^7ri/no/A, n is the total density of a gas, uq is the condensate fraction 
density, and a is the scattering length. Note that the gas parameter no? <^ 1 naturally 
enters the equations. 

First of all, from the Fourier expansion of g{uJn), f{ujn), G{uJn), and -F(w„) we conclude 
that the dissipation in a point- like junction due to the presence of the linear |ci;n| term can 
be associated with the various physical processes. Thus in the case of the g and / contri- 
bution the dissipation can be ascribed to the noncondensate-condensate particle tunneling 
process and exists down to zero temperature [17] . On the other hand, the inspection of the 
expressions for G and F shows that these terms, producing no contribution into dissipation 
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at zero temperature, will be responsible for the explicit T^-behavior of the dissipative effects 
in the junction. The origin of this finite temperature contribution can be attributed to the 
tunneling of thermal phonon-like excitations across the junction. The meaning of the other 
dynamical renormalizations and its temperature behavior will be discussed below. 



IV. JOSEPHSON EQUATION. 

To consider the dynamic behavior of the relative phase </? in the real time, we now follow 
the standard procedure of analytical continuation. Accordingly, the substitution — > —icu 
in the Fourier transform of the Euler-Lagrange equation SSeff/5(f{T) = in imaginary 
time entails the classical equation of motion for the relative phase with the next inverse 
Fourier transformation to the real-time representation. In particular, this means that we 
should replace |a;„| with ~iuj in the above expressions for the Fourier transform of the 
response functions a and (3. The effective action Seff[(f{t)] in the real time, which variation 
5Seff/5(p{t) — yields the real-time equation of motion, can be given by the expression 

, 2 



Seff [(p] = Jdt ^ + 2Io^noinQr cos (p 

a (t-f) cos [(p (t) -(p{t')] 

+P {t - t') cos [(fi (t) + (fi {t')] 



ilSdtdt' < 



(51) 



> . 



In the limit of the slowly varying phase the response functions in the real-time repre- 
sentation can be represented in the form of a functional series. According to Eqs.(48), we 
have 

a{t) = ai5\t) - a25"{t) + a3S"'{t) + ..., 

(52) 

Pit) = PoS{t) - PiS'it) + /32S"{t) - Ps6"'{t) + ■■■. 
Employing variational principle to the effective action Seff[^p{t)] and using decomposition 
(52), we can derive the Josephson equations valid for the slow variations of the phase provided 
the typical time of its evolution is longer than 1/A: 

Gsiip) + (3/2)(/30G'3(^) - ip'Gsiip) + ^G^iip) + {l/2)ip'G',{ip) + ipG{ip) + U'{ip) = 0, (53) 

and in accordance with Eq.(20) 



(f = -Sll{t) ^ III- Il2- 
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Here we have retained the time derivatives of (f{t) to third order corresponding to radiation 
corrections. 

In the following, for the sake of brevity we will consider the case of A; = = A. The 

general expressions for the coefficients for ^ A,, will be considered in the Appendix. 
The potential energy of a junction is given by the well-known relation 



U{f ) = -Ej cos f + {l/2)E2j cos 2ip, 



with the coefficients 



Ej = 2no/o, 



E2J = GoA a - 4i 



TT 



■ 1 1 /27rT\^ 



(54) 



where 



Go = 27/0^ = 167ry — (no/o/A)^ 

A = nou. 



(55) 



In the above relations the condensate density tiq = no(T) depends explicitly on tempera- 
ture according to the well-known expression for a depletion of the condensate fraction of 
the weakly interacting Bose gas [23]. An increase of the temperature leads obviously to 
decreasing the Josephson energy. 

The friction coefficient G{(f) determining the Ohmic dissipation is given by 



G{<p) = Go 



2 1 na-^ f27iT\ . 2 



cos (p+ - 

3 V TT 



sm (p 



(56) 



It is natural that nonzero temperature enhances the energy dissipation of a junction due to 
appearance of thermal excitations in the Bose-condensed gas. 

We may compare Go of the Bose-gas to the analogous conductance of a normal Fermi gas 
with the same density of states: 



Gn = ^TTiyiA), Gn = Sd Go. 



TT 



(57) 



The inverse effective mass of a junction is determined by 



G2{ip) = U ^ - 0:2 - (^2 COS 2(p, 



(58) 
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/32 = 



3G, 



8A 




2 na^ 




(59) 



The renormalization of the effective mass results from the both condensate-noncondensate 
and noncondensate-noncondensate particle tunneling processes. 
The coefficient ^3(92) responsible for the radiation effects reads 



These effects can be associated with emitting a sound from the region of a junction during 
the tunneling of particles across the junction. On the other hand, the radiation effects can 
be treated as a frequency dispersion of the effective friction coefficient. 

We will not enter here in details of the conditions which should be imposed on the 
coefficients of the Josephson equation (53) in order to observe the well-defined Josephson 
effect since this topic is already discussed much in the literature. In essence, the necessary 
condition reduces to the requirement of smallness either quantum zero-point or thermal 
fluctuations for the phase difference across a junction, i.e., mean square value ((A(/?)^) <S 1. 
Involving that {{Aip)"^) ~ T/Ej in the thermal activation region with the crossover at 
To ~ Ej/max{G, ^/G2EJ} to ((A(/9)^) ~ h/max{G, a/G2-Ej} in the quantum fluctuation 
regime at lower temperatures, it is desirable to have sufficiently large Josephson energy Ej 
or, correspondingly, not too large potential barriers. 

Another interesting aspect of such kind of experiment is an investigation of the effects 
beyond the mean- field approximation of a very dilute gas na^ <^ 1, in particular, observation 
of the temperature effects in the dynamics of a junction. As we have seen, the scale of the 
temperature effects should reach the order of yncfi at 2-kT ~ A ^ Tg. This is of much 
interest since the temperature behavior of the Josephson dynamics is closely related to the 
properties of elementary excitations in a condensed Bose-gas. 



tts + /is cos 2lp 




(60) 
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The gas parameter na^ under the conditions typically realized is at most of order of 
10~^, e.g., n ~ lO^^cm"^ and a — 5nm for ^^Rb. In principle, it is possible to approach 
na^ ~ 1, increasing a ^ oo by appropriate tuning of magnetic field with the Feshbach 
resonance. However, too large values of scattering length can facilitate rapid three-body 
recombination. On the other hand, according to recent work [24] such recombination may 
not be inevitable at approaching na^ ~ 10~^. For such values of na^, the effects beyond the 
mean- field approximation becomes well noticeable of about 10%. 

The second interesting aspect is connected with the analog of the voltage-current char- 
acteristic which is an important attribute of a superconducting junction. This kind of 
experiment implies maintenance of the constant bias Sfi for the chemical potentials across a 
junction, corresponding to constant pressure drop 6P = mn6fi and time dependence of the 
relative phase (f{t) — —S/it + ipo. One of possibilities is to use the field of Earth's gravity 
for this purpose. 

As a result, we arrive at the following mean value of the particle current across a Bose 
junction for a sufficiently large period of time as a function of bias S/i 



where 



{G{^)) = -Go 



1 - 



(61) 



1 



ina 



/2nT\' 



and {G^{ip)) — a^- This experiment may give an information on the dissipative properties 
of a junction and also on the nonlinear effects to be of the order of {6n/ A)^. Note that 
nonlinear effects do not contain smallness of (na^)^/^. 



V. CONCLUSION 



To summarize, in this paper we have used a functional integration approach for the model 
of a tunneling Hamiltonian in order to analyze the dynamics of a point-like Josephson junc- 
tion between two weakly non-ideal Bose gases. The effective action and response functions 
which describe completely the dynamics of a junction are found. Using the low frequency 
decomposition of the response functions, the quasiclassical Josephson equation which the 
time evolution of the phase difference (p across the junction obeys is obtained to the terms 
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of third order in time derivatives. The corresponding kinetic coefficients are calculated ana- 
lytically, involving the finite-temperature corrections. The temperature effect on the kinetic 
coefficients demonstrates the T^-behavior. 

Like a junction of the planar geometry [15, 17, 22], the dynamics of a point-like junction 
concerned here has an Ohmic dissipative type due to the gapless character of excitations in 
a Bose-condensed gas. Note only that the scale of the dissipative effect in the latter case is 
significantly less. 

The behavior of the dissipative particle current on the temperature and the difference 
of the chemical potentials in Eq. (61) can be used to investigate the role of quasiparticle 
excitations in the dynamics of a Bose junction. Since the dissipation is closely related to 
the tunneling process of condensate particles and depends on the structure of a junction, 
the effect of dissipation in the junction dynamics can be reduced provided the structure of a 
junction prevents condensate particles from tunneling across the junction. We believe this 
question deserves further study. 
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VII. APPENDIX 

Here we present the general case of A; 7^ A^. To derive expressions for the coefficients in 
the equation for the phase (53), we first find the Taylor expansion of the response functions 
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a and (5. Finally, for the coefficients in the functional series (48) we arrive at the expressions: 



as 



2 1 ir- 3 
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•A 1/4 
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^^3 — TeA^A^ '^^r ~^ 15 

Here we use the following notations 
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and (j) is introduced in (34). 

Then, for the potential energy of a junction we have 



Ej = 2Ioy^noinor, 



3 3 \ 1/4 

nia'^nrcij. ^ 
^2 



In 



4a;e 
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Go is given by 



Go = 27/0 = 167r 



3 3\ 1/4 / / 



^ V 

Ai,r = UQirULr- 



(62) 



(63) 



(64) 



(65) 



In the above relations the condensate densities no — no{T) depend on temperature according 
to the well-known expression for the depletion of the condensate fraction of the weakly 

interacting Bose gas [23]. 

The friction coefficient G{(p) determining the Ohmic dissipation is represented by 

2 



G{ip) = Go 



bil^ cos^ ip + 



1 / 3 3\ 1/4 

1 / niaiUra'^ \ 



27rr 



sin^ ip 
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It is natural that nonzero temperature enhances the energy dissipation of a junction due to 
appearance of thermal excitations in the Bose-condensed gas. 

Comparing Go for the Bose-gases with the conductance of normal Fermi gases with the 
same densities of states, we have 



(66) 



The inverse effective mass of a junction is determined by 



^0 J 1,(3) 



3 3X1/4 



P2 



3Gn 



d3) 



2 
3 



/ o 1/1 



0i2 - P2 COS 2(/7, (67) 
2" 




(68) 



The renormalization of the effective mass results from the both condensate-noncondensate 
and noncondensate-noncondensate particle tunneling processes. 
The coefficient Gs{ip) responsible for the radiation effects reads 



^0 I , (4) 



3 3\ 

niafrira^ \ 
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1 - 



16A/Ar VVAIA; 
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16AiAr 
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(69) 



The expression derived for the dynamical coefficients in equation (53) governing the phase 
difference across a point-like junction allows us to describe the low frequency dynamics in 
the asymmetric case of Bose gases with the different order parameters. 
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